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Abstract: We determine the S-matrix that describes scattering of arbitrary bound states 
in the light-cone string theory in AdSs x S 5 . The corresponding construction relies on the 
Yangian symmetry and the superspace formalism for the bound state representations. The 
basic analytic structure supporting the S-matrix entries turns out to be the hyper geometric 
function 4F3. We show that for particular bound state numbers it reproduces all the 
scattering matrices previously obtained in the literature. Our findings should be relevant 
for the TBA and Liischer approaches to the finite-size spectral problem. They also shed 
some light on the construction of the universal R-matrix for the centrally-extended psu(2|2) 
superalgebra. 
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1. Introduction 

The S-matrix approach has been recently recognized as an indispensable tool to study the 
spectral problem of both the AdSs x S 5 superstring and the dual gauge theory [1-10]. 
In the uniform light-cone gauge [11, 12] the corresponding string sigma model describes a 
two-dimensional massive quantum field theory of eight bosons and eight fermions. The 
light-cone momentum P + , which is a gauge fixing parameter, plays the role of the string 
length. In the limit when P + tends to infinity, the sigma model is defined on a plane, which 
calls for the application of scattering theory. 

For a scattering process, the statement of integrability implies absence of particle pro- 
duction and factorization of multi-particle scattering into a sequence of two-body events. 
Since the sigma model does not have two-dimensional Lorentz invariance, the two-body 
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S-matrix for scattering of fundamental particles (giant magnons [13]) has a rather intricate 
structure. It is invariant w.r.t. the superalgebra psu(2|2) © psu(2|2) enhanced by three 
central charges depending on particle momenta [6]. The latter algebra is a symmetry of 
the light-cone gauge fixed Hamiltonian in the off-shell theory where the level-matching 
condition, i.e. the requirement of vanishing of the total world-sheet momentum, is sus- 
pended [14]. The matrix structure of the two-body S-matrix appears to be uniquely fixed 
by the centrally-extended psu(2|2) © psu(2|2) algebra, the Yang-Baxter equation and the 
generalized physical unitarity condition [6, 15, 16], while the overall scaling factor is severely 
constrained by crossing symmetry [17]. 

In the limit of infinite light-cone momentum, in addition to the fundamental particles, 
the spectrum of the string sigma model contains an infinite tower of bound states [18]. 
More explicitly, ^-particle bound states comprise into the tensor product of two 4£-dim 
atypical totally symmetric multiplets of the centrally-extended psu(2|2) algebra [19]. 

Concerning the finite P+ spectrum, any power- like 1/P+ corrections to the energy 
of a multi-particle state can be obtained by means of the Bethe-Yang equations [5, 20, 21] 
representing the quantization condition for the particle momenta. A complete handle on the 
string asymptotic spectrum and the associated Bethe-Yang equations requires, in principle, 
the knowledge of the psu(2|2)-invariant S-matrices § £l ^ 2 which describe the scattering of 
i\- and ^-particle bound states. 

In addition to the power-like corrections to the spectrum, there are also exponentially 
small corrections which are not captured by the Bethe-Yang equations. The leading expo- 
nential corrections [22-24] to the dispersion relation for the fundamental particles and the 
bound states can be derived by applying the perturbative Liischer's approach [25], which is 
also based on the knowledge of the world-sheet S-matrix [26-30]. It appears that Liischer's 
approach could also be applied to find perturbative scaling dimensions of gauge theory 
operators up to the first order where the Bethe-Yang description breaks down [31, 32] (see 
also [33]). The energy of the two-particle state corresponding to the Konishi operator re- 
ceives contribution from the whole tower of ^-particle bound states which travel around a 
cylinder of finite circumference P+. The corresponding computation thus exploits the scat- 
tering matrices of fundamental particles with bound states, S u . A more general treatment 
of the scaling dimension of a gauge theory operator corresponding to an on-shell bound 
state particle will thus require the knowledge of a generic scattering matrix § £l ^ 2 . 

To find the exact, i.e. non-perturbative string theory spectrum, one could try to adapt 
the Thermodynamic Bethe Ansatz (TBA) approach, originally developed for relativistic 
integrable models [34]. To this end, one has to find the scattering matrices for bound 
states of the accompanying mirror theory [16]. The bound states of this theory are "mirror 
reflections" of those of the original string model. Furthermore, the scattering matrix of 
the mirror particles is obtained from §^ 2 by double Wick rotation and the transfer to 
the anti-symmetric representations. Although the Bethe-Yang equations for bound states 
of the original and the mirror model can be obtained by fusing the Bethe equations for 
the corresponding fundamental particles 1 , the knowledge of §^ 2 might be relevant for 
an alternative approach based on functional relations between eigenvalues of the transfer 
1 We refer to [35] for a derivation of the corresponding Bethe equations from Yangian symmetry. See 
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matrices. In particular, it would be interesting to check various conjectures about such 
eigenvalues in higher rank representations. For an interesting recent development in this 
direction we refer to [39-42]. 

The aim of this paper is to complete the program of [43] by finding the S-matrix 
g^ife w hich corresponds to scattering of bound states with arbitrary bound state numbers. 
This will be done by using the Yangian symmetry [44] in conjunction with the superspace 
formalism [43]. 

According to [43], the ^-particle bound state representation % of the centrally extended 
su(2|2) algebra can be realized on the space of homogeneous (super) symmetric polynomials 
of degree t depending on two bosonic and two fermionic variables, w a and 6 a , respectively. 
Thus, the representation space is spanned by an irreducible short superfield $>e(w, 9). In 
this realization the algebra generators are represented by differential operators linear in JJ 
in w a and 6 a with coefficients depending on the representation parameters (the particle 
momenta). The S-matrix §^ 2 acts on the product of two superfields <J>^ (w , 6)<&£ 2 {v , d) as 
a differential operator of degree i\ + £2- We require this operator to obey the following 
intertwining property 

S^ 2 • J12 = J21 • , 

where J12 and J21 are the generators of the centrally extended psu(2|2) in the corresponding 
two-particle representation. For su(2) subalgebras of psu(2|2) this condition literally means 
the invariance of the S-matrix, while for the supersymmetry generators it involves the 
braiding (non-local) factors [15, 45, 46] to be discussed later. In [43] the lower-dimensional 
examples S 11 , § 12 and 8 22 have been found by solving the above invariance condition 
together with the Yang-Baxter equation. The fundamental S-matrix S 11 has been shown 
[44] to commute with the Yangian for the centrally extended psu(2|2) algebra. It appears 
that the new examples § 12 and § 22 also respect Yangian symmetry which provides an 
alternative to solving the Yang-Baxter equation [47] . Denoting by J a Yangian generator 
acting in the finite-dimensional evaluation representation corresponding to a bound state, 
we thus require 

§^ 2 • J12 = J21 • §^ 2 , 

for any 8^ 2 . The Yangian carries a natural Hopf algebra structure [44], so that J12 and 
J21 can be regarded as the coproducts A (J) and A op (J), respectively, evaluated on bound 
state representations. Here A op stands for the opposite coproduct. 

Our construction of §^ 2 will involve three two-particle bases: the first is the standard 
one corresponding to the product ^ (w, 0)&e 2 (v , $), the second is given by certain prod- 
ucts A(JI ai ) • • • A(J°fc) and the third by A°P(S ai ) ■ ■ ■ A°P0 ak ). Denoting by A and A°p the 
transition matrices from the second and the third basis to the standard one, we will show 
that the matrix representation of the S-operator in the standard basis is nothing else but 

S* 1 * 2 = A op ■ A -1 . 

the references in [36] for an earlier literature on Yangians in AdS/CFT. Recent progress can be found in 
[37, 38]. 
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Finding A and A op , we will thus establish the matrix form of §^ 2 . As we will see, the 
basic analytic structure supporting the expressions for the corresponding matrix elements 
is given by the hyper geometric function 4F3. We point out that our construction of § 
involves linear algebra only and it can be easily implemented in the Mathematica program. 
It is interesting to point out that the above matrix factorization of the S-matrix reminds 
of the Borel decomposition of the double of the Yangian. 

Taking the tensor product of two psu(2|2)-invariant S-matrices, one obtains psu(2|2) © 
psu(2|2)-invariant world-sheet S-matrices which describe the scattering of bound states in 
the light-cone string theory on AdSs x S 5 . Of course, the resulting S-matrix should be 
multiplied with the proper overall scalar factor obeying crossing symmetry. This factor 
has been already obtained in [43, 48, 49]. 

Having found the general bound state scattering matrix E> £l£2 , we then verify that it 
reproduces all the previously obtained special cases and satisfies the necessary physicality 
requirements. In particular, we investigated its analytic structure and confirmed that it 
only exhibits the expected physical pole corresponding to the formation of a bound state 
of rank t\ + £2 ■ This completes the bound state S-matrix program. 

It is worth stressing that, besides the uses in the TBA or Liischer's approaches , which 
are physically the most relevant ones, there is also a mathematical interest in deriving 
the universal R/S-matrix for the centrally-extended psu(2|2) algebra. The solution to this 
problem has been so far elusive, in spite of the progress nevertheless achieved [50-55]. The 
present work provides an important step towards finding the universal R-matrix. In fact, 
our formulae present the explicit R-matrix entries for arbitrary evaluation representations 
of the relevant Yangian, namely, all (finite-dimensional) bound-state representations. On 
one hand, this exhausts a large class of representations. On the other hand, it may now 
become easier to establish a universal form of the R-matrix which could reproduce all these 
expressions, as started in [36]. 

The paper is organized as follows. In the next section we introduce our notations and 
formulate our main result - the S-matrix §^ 2 of arbitrary bound states in the standard 
basis. The rest of the paper contains details of the derivation and verifications of various 
properties. Discussion of certain analytic aspects is relegated to the appendices. 



2. Kinematical Structure of the S-Matrix 

In this section, we will discuss the kinematical structure of the S-matrix. In particular, we 
will use su(2|2) invariance to show that the S-matrix is of block diagonal form. 

2.1 Centrally extended su(2|2) 

We will first discuss centrally extended su(2|2). This algebra has bosonic generators R,L, 
supersymmetry generators Q, G and central charges H, C, C'. The non-trivial commutation 
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relations between the generators are given by 



[L a b , J c ] — 5 b c $ a — \8\j c , [IRaf, Xy] — SyS a ~ ^aS-y, 

{Q a a , Qfl 6 } = e^C, {G a a , G/} = e^e a6 Ct, 1 ' j 

{Q„, Gf} = ^Rof + ^ V + ±«S«<^ EL 

The eigenvalues of the central charges are denoted by H, C, . The charge M is Hermitian 
and the charges C, and the generators Q, G are conjugate to each other. 

For computational purposes, it proves worthwhile to consider representations of the 
algebra in the superspace formalism. Consider the vector space of analytic functions of two 
bosonic variables w\ t 2 and two fermionic variables #3^. Since we are dealing with analytic 
functions we can expand any such function $(w,0): 

00 

*(w,9) = 0), 
£=0 

$ € = cj) a ^ Wai ... Wat + <t> ai ~ a '- ia w ai . . . w ai _ x e a + 

f- fl » a V... %2 fl^. (2.2) 

The representation that describes ^-particle bound states is 41 dimensional. It is realized 
on a graded vector space with basis \e aim .. ag ), |e 01 ... a£ _ lQ ), |e ai ... 0£ _ 2Q( g), where ai are bosonic 
indices and a, j3 are fermionic indices, and each of the basis vectors is totally symmetric 
in the bosonic indices and anti-symmetric in the fermionic indices [18,19,43]. In terms 
of the above analytic functions, the basis vectors of the totally symmetric representation 
can evidently be identified as \e ai ... ai ) <-»■ w ai . . . w av \e ai ... at _ ia ) w ai . . .w ai _ x Q a and 
\ e a 1 ...a e _ 1 a8) ^ w ai . . . w ae _ 2 9 a 9p, respectively. In other words, we find the atypical totally 
symmetric representation describing ^-particle bound states when we restrict to terms $£. 

In this representation the algebra generators can be written in differential operator 
form as 



jf 6 d ltd,,, d up ft n d lA^fl 9 

— w aQ Wb 2°a w cQ Wc , —VoL dd(j 2 0a0 l 89-,' 

Qa" = + be ab e a pw b 4-, G a « = dw a ^ + ee ab e«%^ 



(2.3) 



and the central charges are 



H 



ab (w a ^ + 6 a JL.) , Ct = ed + e a JL.) , 

(ad + bc) (w a ^- + e a ^y 



(2.4) 



To form a representation, the parameters a, b, c, d must satisfy the condition ad — be = 1. 
The central charges become £ dependent: 

H = £{ad + be), C = lab, C f = led. (2.5) 
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The parameters a, b, c, d can be expressed in terms of the particle momentum p and the 
coupling g: 



„ _ / 9 n h _ / 9 K ( 

r — — Hi 77 rl — /lExl. ( -\ — x ~ 

c — A/ 2££3F' ~~ V 2i ir > V ^ 



(2.6) 



where the parameters x 1 * 1 satisfy 



x+ + -1 - x- - — = — , — = e i 'P (2.7) 
x + x g x 

and the parameters r] are given by 

rj = e l ^rj(p), t](p) = e l ^\Jix~ - ix+ , ( = e 2 ^. (2.8) 

The fundamental representation, which is used in the derivation of the S-matrix scattering 
two fundamental multiplets [6, 15], is obtained by taking 1=1. 

The bound state S-matrices should respect the su(2|2) symmetry, by requiring invari- 
ance under the coproducts of the generators 

S A(J A ) = A op (J A ) S, (2.9) 

where A op = PA, with P the graded permutation 2 . 

2.2 Invariant subspaces 

Consider two bound states with bound state numbers £1,^2 respectively. The tensor prod- 
uct of the corresponding bound state representations in superspace is given by: 

*(u;,0)S(t;,0), (2.10) 

where w, 6 denote the superspace variables of the first particle and v, 1? describe the repre- 
sentation of the second particle. 

The S-matrix acts on this tensor space and should, according to (|2.S|), commute with 
A(L)} and A(M)3. Prom this, it is easily deduced that the numbers 

K 11 = #0 3 + #0a + ##3 + #^4 + 2#w 2 + 2#v 2 , 
K m = #0 3 + #^3 + #w 2 + #v 2 (2.11) 

are conserved. The variables W2,V2 can be interpreted as being a combined state of two 
fermions of different type [6]. Hence, the number K 11 corresponds to the total number of 
fermions, and the number K 111 counts the number of fermions of one species, say, of type 
3. The fact that these numbers are conserved allows us to define invariant subspaces, for 
each of which we will derive the corresponding S-matrix. 



2 We remind that, in the non-local formalism of [15], one needs to explicitly permute the two represen- 
tations when acting with P. 
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Let us write out the tensor product more explicitly. Since we are considering bound 
states with bound state number £±,£2 we restrict to 



h-k-X k-X 

w 2 

l-Xl-U 



"2 + W""X ) 



(2.12) 



The ranges over which the labels k, I are allowed to vary can be straightforwardly read off 
for each term. By multiplying everything out, we reproduce all the basis vectors. One can 
compute the quantum numbers K ,K for any of these basis vectors. The results are 
listed in Table |]. When we take a closer look at the result, we see that when we order the 



Space 1 



Space 2 



K 



in 



N 



Case 



h-k-Xk 



h-i-Xi 



3 w{ 



h-k- 



1 fr 



2(k + l)+2 



h-l-x 1 

v 2 



-k-X k 
Wo 



-L.l 



w 



1 



e 3 w\ 
e 3 e 4 w{ 



W 2 

-k-Xk 



W 2 

-k-X k-X 
Wo 



$A 2 



h-l-X n l 



h-i-x l-X 

1 u 2 

-i-x„.i 



-k- l „,k 



Wo 



-K,i 



w 



1 

9 A w\ 



Wo 



i-k-X k 
W$ 



a a £\—k—X k- 



$±v[ 2 



<a-i-l.,J 



2(fc + Z 

2(k + l 
2(k + l 
2(k + l 
2(k + l 

2{k + l 
2(k + l 
2{k + l 
2{k + l 







v e r l v i 2 


w[^ k w k 






e 3 e^~ 


k - l w k - 1 


V2 2 ~ l V l 2 


e 3 e A w[^- 


k ~ l w k ~ l 




e 3 w"r k - 


l w k 




e A w l r k - 


l w k 





J2-1-X l-X 

°2 

v 2 
v 2 



2(k + l) 
2{k + l) 
2{k + l) 
2{k + l) 
2(k + l + l) 
2{k + l + l) 



+ 2 

+ 1 
+ 1 
+ 1 
+ 1 

+ 1 
+ 1 
+ 1 
+ 1 



k + l + 2 

k + l 

k + l + 1 
k + l + 1 
k + l + 1 
k + l + 1 

k + l 
k + l 
k + l 
k + l 

k + l 
k + l 
k + l 
k + l 
k + l + 1 
k + l + 1 



k + 

k + 

k + 
k + 
k + 
k + 

k + 
k + 
k + 
k + 

k + 
k + 
k + 
k + 
k + 
k + 



+ 1 
+ 1 



la 

lb 

Ha 
Ha 
Ha 
Ha 

lib 
lib 
lib 
lib 



Table 1: The 1%1\£2 vectors from the tensor product and their su(2) x su(2) quantum numbers. 

states by the quantum numbers K ll ,K 111 , there are exactly five different types of states: 

Case la: K u = 2N + 2, K m = N + 2, 

Case lb: K n = 2N + 2, K m = N, 

Case Ha: K u = 2N + l,K m = N + 1, 

Case lib: K 11 = 2N + 1, K m = N, 

Case III: K n = 2N, K m = N, 
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for some integer N. Each of these states have different quantum numbers K ,K , hence 
the states belonging to each of these cases form a subspace which is left invariant by the 
S-matrix. 

Clearly, vectors from Case la and Case lb only differ by the exchange 3 4, which 
is easily realized in terms of the (fermionic) su(2) symmetry generators of type R. Hence, 
the subspaces spanned by the two types of states are isomorphic, and scatter via the same 
S-matrix. An analogous relationship connects Case Ha and lib. Thus, there are only three 
non-equivalent cases: 

Case I: K u = 2N + 2, K m = N + 2, 

Case II: K n = 2N + 1, K m = N + 1, 

Case III: K u = 2N, K m = N. 

For fixed N (i.e. for fixed K ll ,K 111 ) we denote the vector spaces, spanned by vectors from 
each of the inequivalent cases, by V , V 11 , V respectively. 

2.3 Basis and relations 

Later on we will introduce different bases for the different cases, but in this section we will 
discuss the basis as obtained by multiplying out the tensor product as seen from Table [j]. 
We will call this basis the standard basis. 

Case I, K n = 2N + 2, K m = N + 2. 

For fixed N, the vector space of states V 1 is N + 1-dimensional. The standard basis for 
this vector space is 

\k, I) 1 = Ozw[^- l w k z ^svt^M , (2.13) 

Spacel Space2 

for all k + I = N. These indeed give N + 1 different vectors. 
Case II, K u = 2N + 1, K m = N + 1. 

For fixed N, the dimension of this vector space is 4iV + 2. The standard basis is 

\k,i)Y = e-^-^wl vf- l v{, 

= M 2 ~'~H (2-14) 



\k,i)^ = e.e^-^wt 1 #svi 2 - l - 1J 



2' 



where k + I = N . As a lighter notation, we will from now on, with no risk of confusion, 
omit indicating "Space 1" and "Space 2" under the curly brackets. The ranges of k, I are 
clear from the explicit expressions and it is easily seen that we get 4iV + 2 states. 
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Case III: K u = 2N, K m = N 

For fixed N = k + I, the dimension of this vector space is 6N. The standard basis is 



I*, Of 1 

\k,l)f 



™1 w 2 y v l v 2> 



u 2 i 



IM> 



in 

3 



IM> 4 n 




-Z-l Z-l 
u 2 ' 



(2.15) 



in 



\k,l) 5 
IM>F 



k w k - x 



v 2 



Note that our numbering slightly differs from the one used in Table [T], in the sense that 
\k, are rescaled in such a way that they also have N = k + I, instead of A; + I + 1 as in 
Tabled 

It is convenient to supply all these spaces with the usual inner product 3 



f(k,l\m,n)t 



(2.16) 



We also introduce the vector spaces = span{|M)f}, for A = I, II, III. Note that for 
M/Owe have dim = 1, dim V^j = 4, dim = 3, dim = 2 etc. 

The different cases are not unrelated. One can use the (opposite) coproducts of the 
symmetry generators to move between the different subspaces. In particular, the cases are 
distinguished by their quantum numbers K 11 , K 111 . Acting with supersymmetry generators 
will change these numbers. Hence, these generators provide maps between the cases. How 
this works is schematically depicted in Figure |]. 

These relations between the different cases will play an important role in the derivation 
of the full S-matrix. In the next section we will introduce two different sets of bases which 
allow for a natural interpretation of the S-matrix. These bases will make use of the full 
Yangian symmetry rather than just the su(2|2) as we did in this section. In this framework 
we can solve Case I. Then we employ the different arrows in Figure |l] (and their Yangian 
counterparts) to relate the different S-matrices to the Case I S-matrix. 

Summarizing, we find that the S-matrix is of block-diagonal form 



/ X 



W\ 







\ 







<3f 



% ) 



(2.17) 



3 For the purpose of the following derivation, one can convince oneself that orthogonality of these vectors 
is actually sufficient. 
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Figure 1: Schematic representation of the relations between the invariant subspaces. The opposite 
coproducts also respect the above diagram, as well as all their Yangian counterparts. 



The outer blocks scatter states from V 1 



SC : V L — ► V L (2.18) 

k+l 

IM) 1 — E &m l \m,k + l- m)\ (2.19) 

m=0 

where &m will be given by ( 4.1 1| ) . The blocks & describe the scattering of states from 
V 11 

W : V 11 — > V u (2.20) 

k+l 4 
m=0 j=l 

These S-matrix elements are given in ( |5.18p . Finally, the middle block deals with the third 



case 



& : V m — > V UL (2.22) 

k+l 6 

i fc > of -EE "'■ k+l ~ ™>f. ( 2 - 23 ) 

m=0 j=l 

with from doil) . 

We recall that the full AdSs x S 5 string bound state S-matrix is then obtained by 
taking two copies and multiplying each one of them with the phase factor [43] 

<Sb(Pi,J»2) = f^) f^rl cr(xi,x 2 )x 

4-1 

x VG(^ 2 - ii)G(£ 2 + h) J] + 2 ^)> (2-24) 

9=1 
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where, in our conventions, 



G(Q) _ *Z*±$ (2.25) 



gives the standard pole at bound-state number Q (see Section 4L2 later on for a discussion 
on the physical poles of the S-matrix) . Here u is given by 

u = IT- \ x+ + -X + x ~ + — J • ( 2 - 26 ) 

3. Yangian Symmetry and Coproducts 

So far we have only used su(2|2) symmetry to study the bound state S-matrix. This, 
however, is not enough to fix the S-matrix (up to an overall phase). In particular, it was 
found that one needs to impose the Yang-Baxter equation by hand to attain this [43]. An 
alternative to this method was shown to come from Yangian symmetry [47]. 

The Yangian of su(2|2) can be supplied with a Hopf-algebra structure. In this language 
the invariance of the S-matrix can be formulated as 

§ A(J A ) = A op (J A ) §, 

§ A(P) = A op (J A ) §, (3.1) 

where J stands for any generator of the Yangian. For explicit formulae and more details 
we refer to Appendix [A|. All this seems to indicate that the full Yangian of su(2|2) should 
be viewed as the underlying symmetry algebra of scattering processes. Indeed, as we will 
see later on, we are able to construct any bound state S-matrix from this algebra. 

3.1 (Opposite) coproduct basis 

Let us turn back to the invariant subspaces. In this section, we define different bases for 
each case apart from the standard basis, which is commonly used in the literature. We 
call them the coproduct basis and the opposite coproduct basis. The basis transformation 
between the coproduct (opposite coproduct) basis and the standard one will be denoted 
by A (A op , respectively). 

These bases will be constructed by using Yangian generators to create states out of a 
chosen vacuum. This is similar to [35] where it was used to study the Bethe Ansatz. We 
define our vacuum to be 

|0> = wM 2 , (3.2) 

just as it is used in the coordinate Bethe Ansatz. We normalize our S-matrix in such a 
way that S| 0) = |0). The (opposite) coproduct basis will consist of states created by the 
(opposite) coproducts of various symmetry generators acting on this vacuum. Clearly, the 
S-matrix has a natural interpretation in these bases, and can be formulated in terms of A 



and A op , as will be explained in section 3^. We will now list the explicit formulae for the 
different bases. 
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Case I, K 11 = 2N + 2, K m = N + 2. 

The coproduct basis is given by 

N q 

A(Qi)A(G|) J] A(Li)fjA(Li)|0), q = 0, 1, . . . N, (3.3) 

i=q+l j = l 

and the opposite coproduct basis is given by 

N k 

A<*(Qi)A<*(G!) I] A^(Li)nA op (Li)|0), fc = 0, 1, . . .iV. (3.4) 

i=k+i j=i 

Each of these two bases is indeed composed of N+l different vectors. By explicitly working 
out the coproducts one can see that these vectors form a basis for Case I. One could also 
consider an alternative choice, like for instance 

N k 

A(Q 3 )A(Q 3 ) [] A(I4);QA(£4)|0), (3.5) 

i=k+l j=l 



but these vectors are readily seen to be proportional to (3.3). 



It is also straightforwardly seen why Q3,3j ) actually describes Case I from the point 
of view of the quantum numbers K , K . The operators AL^AL^ create a boson of 
type 2 out of the vacuum and the supersymmetry generators AQ 3 , AQ 3 create a fermion 
of type 3. Hence we find that K u = 2#L^ + 2#AL£ + #AQ^ + #AQ^ and K m = 

+ #AL| + #AQ^ + #AQ\. This indeed coincides with K u = 2N + 2, K m = N + 2. 

Case II, K u = 2N + 1, K m = N + l. 

The coproduct basis is given by 

N q 

A(Q 3 ) [] A^nACL^O), 

i= q+ X j=i 

N q 

A(Ql) J] A^HA^IO), 

i=q+l j = l 

N-l q 

A(Qi)A(Qi)A(Qi) J] ACL^IJAOL^IO), 

i=q+l j=l 
N-l q 

A(Q 3 )A(Q 3 )A(Qi) H A(Li)fjA(Li)|0), (3.6) 

i=q+l j=l 

and similar expressions hold for the opposite coproduct basis. One can again compute 
K ,K for these states and see explicitly that they describe Case II. 
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Case III, K n = 2N, K m = N. 

The coproduct basis is 

N q 

n A(Li)nA(£4)io>, 

i=q+l j=l 

N-l q 

A(Qi)A(Q|) J] A(l4)nA(q)|0), 

i=g+l 3=1 

AT-1 q 

A(Qi)A(Qi) [J A(Li)nA(Li)|0), 

«=<?+i i=i 
AT-l <? 

A(Qi)A(Qi) [J A(Li)nA(Li)|0), 

i=g+l j = l 

AT-1 q 

A(Qi)A(Ql) [] A(l4)IjA(]Li)|0), 

N-2 q 

A(Qi)A(Qi)A(Q|)A(Qi) [] A(L^) A(L^)|0). 

i=9+l i=i 

These are readily seen to be 6^ states and their quantum numbers are of the form K = 
2N, K m = N. 

The Yangian generators also provide maps between the different cases. In particular, 
one finds that Figure [j] also holds for Yangian generators. One important thing to notice 
is the following. Even though, for example, AQ| maps Case II onto Case I, this does not 
automatically give a straightforward map between the vector spaces Vu). For instance, one 
has 

AQ3 : Vjjij — > Vlj_ x © Vl_ x>l , (3.7) 
ml : Vl] — » ^Vl^-l © Ki © VU,1+1- (3-8) 

These relations between the different cases will be used later on. 

3.2 S-matrix in coproduct basis 

The fact that the coproduct basis is well suited for computing the S-matrix can be seen 
from (p.l|). One sees that the S-matrix directly maps the coproduct basis onto the opposite 
coproduct basis. In particular, since we normalize the S-matrix in such a way that S|0) = 
|0), we see that the S-matrix, when written as a map between these two bases, is just the 
identity matrix. 

In other words, one can now get the general formula for the S-matrix in the standard 
basis just by applying the appropriate basis transformations. Let us denote the S-matrix 
written in the standard basis as S. One then finds: 

S = A op A~ 1 . (3.9) 
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Note that the explicit matrices A and A op respectively, just consist of the coproduct vectors 
written in the standard basis. The above discussion can be summarized in the following 
commutative diagram: 



{coproduct basis} > {opposite coproduct basis} 



A 



(3.10) 



{standard basis} ► {standard basis}. 



The computationally hard part is finding the explicit inverse of A. For any explicit case at 
hand this can be done by simple linear algebra, but the expressions become rather involved. 
However, we will be able to carry out this procedure in full generality for the S-matrix of 
Case I, and use this result to find the S-matrix for all the other states. To illustrate the 
above discussion, we conclude this section with an explicit example. 

One of the easiest examples is that of K 11 = 1, which can be considered as a warm- 
up for Case II. The vector space consisting of these states is only two dimensional. The 
standard basis vectors are 

M 1_ M a . ti^gt^- 1 . ( 3 - n ) 
The coproduct basis is given by 

A(Qi)|0), A(Qi)|0). (3.12) 
Explicitly, when written down in terms of the standard basis, this gives: 



v(Qi)|0) = MiM^M 2 + a 2 £ 2 w{^ 3 v{ 



e 2 -i 



A(Q 3 )|0) = ^{^(aida-^Xa-iam}^,!^ + (3 ^ 

, h{(ai(hc2 - a 2 di)£ 1 - ia 2 u 2 } tl „ e 2 -i 
+ ^ tOiWaV , 

or, more conveniently written in matrix form, 

( a £ £i{(a2(aid2-b2Ci)e2-iaiui} \ 
fl2 ^ 2 l2{C,ai(-a2di+b 1 c 2 )l 1 -ia2U2} J ■ ( 3 - 14 ) 

One obtains a similar expression for the basis transformation concerning A° p . 

It is now an easy exercise to compute § for this case via ( |3.9[) . Doing the algebra gives 
the following result: 



\ \fhr]{v2) xf-x 2 xf-x 2 / 

which indeed coincides with the known S-matrices [6, 15,31,43], and also with the coeffi- 
cients recently found in [35]. 



.E2. X+-X+ jfy y/h V (p 2 ) x+- Xl 

V2 
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4. Complete Solution of Case I 



In this section, we will discuss the S-matrix of Case I . We will use Yangian symmetry to 
derive its explicit form. This S-matrix proves to be the building block out of which the 
S-matrices for both Case II and Case III can be constructed 4 . Because of this, we will study 
it in detail before moving on to the other cases. We will analyze its pole structure and also 
compare this S-matrix in the semi-classical limit to the classical r-matrix proposed in [52] . 

4.1 Explicit solution 

Our starting point is the coproduct basis (|3.3|), It is convenient to reorder the products in 
the following way: 

N q I 

I] A(Ljjj) A(L 2 ) \ A(Qi)A(G 4 )|0). (4.1) 

J=q+1 j=l J 

The action of the susy generators on the vacuum is of the form 

A(Qi)A(G 4 )|0) = (a 2Cl - a lC2 )^ 2 |0, 0)\ (4.2) 

with a similar expression for the opposite version. As a matter of fact, from this one can 
straightforwardly read off the action of the S-matrix on |0, 0) 1 : 

, _ 8A(Qj)A«a|0> 
(a 2 ci - aic 2 )«i«2 
A°P(Ql)A op (Gj)§\0) 
(a 2 ci - aic 2 )£ih 



- a 3 c 4 
a 2 ci — aic 2 



CO) 1 . (4.3) 



In other words, the S-matrix multiplies |0, 0) 1 by a scalar. We will denote this scalar by D. 
In terms of x^, it is given by 

D _ a 4 c 3 - Q3C4 _ £1 - e 2 ^ 44 ^ 
a 2 ci — aic 2 xf — X2 e i? z 

It is readily seen that this coefficient is indeed consistent with the action of the S-matrices 
previously found in [6, 15, 31, 43]. 

One can now use the generators L 2 ,L 2 to construct a generic Case I state \k,l} 1 from 
|0, 0) , for arbitrary k, I. This can be seen by considering the following identities: 

(L 2 l)(5u + A(h\))\k, I) 1 = {A(L 2 ) - u 2 A{lL 1 2 ) + A(L 2 ) o (h\ ® 1)} \k, I) 1 , (4.5) 
(1 ® L^)(5u + A(L\))\k, I) 1 = - |a(i4) - u a A(L 2 ) - A(L^) 0(1® L})| \k, l}\ 



4 Our procedure will somehow be reminiscent of employing highest weight states of Yangians. 
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where 



5u = U\ — U2- 



(4.6) 



Since A(L})|/c, I) 1 = ^ 1+ ^ 2 ~^ fc+ ^ +1 ) it is obvious that the left hand side is proportional 

to \k + 1,1) and \k,l + 1) respectively. By applying these operators inductively to |0,0), 
one finds 



k+l 



\k,l) 



m=l n=l g=l ^ 

[(hi ® l)(6u + A(h\))} k [(1 ® hl)(Su + A(L}))] ' |0, 0). 



Then, by (4.5) 
\k,l) l = 



(4.7) 
(4.8) 





\A(Ll) + ^- 2u2 2 - 2i+1 A(L|)l 








Ut=i(h - m) ni=i(^2 - n) U k q t[ (Su + ^- q) 





ICO) 1 



This exactly tells us how to write a state in the standard basis as a combination of coprod- 
ucts. In other words, this explicitly indicates how to construct A , It is now straight- 
forward to obtain the action of the S-matrix on Case I states from this. The symmetry 
properties of the S-matrix, together with (|4.3| ), now imply 

SIM) 1 = 



D- 



nti 






|"2ui+«a-l-!y A°P(L|) A°P(Ll)l 









(4.9) 
I0.0) 1 . 



The right hand side can be computed straightforwardly. One finds that SjA;,/) 1 is of the 
form 

k+l 



with 



k.l 



D 



n=0 



n 



I 



rt i (h -p)uU (£2 - P ) nS (su + « - P ) 



(4.10) 



(4.11) 



n—m 



k I / L \ / 7 \ m l—n k—m 

*e (*_ m )(J m )n«f n?n«^n« 

m=0 I x 7 x 7 p=l p=l— m p=l 



fc-H — m — p+2 
2 -»■ 2 



The coefficients are given by 



y in 



<Jit± 



<5u± 



2 

^1 + ^2 



m + 1, 
m + 1, 



(4.12) 



£i + l- 2i, 
£2 + 1- 2i. 
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It is worthwhile to note that in the special case I = (and similarly for k = 0) this 
expression reduces considerably. For later use, we can write it in the following way: 



k ~ n \nj TT* iu+^-i, 



In all of the above expressions it is understood that products are set to 1 whenever they 
run over negative integers, i.e. Yl a = 1 if 6 < a, an d the binomial (^) is taken to be zero if 
y > x and if y < 0. 

We can notice how the formula we have found bears a rational dependence on the dif- 
ference of the spectral parameters, as typical of Yangian universal R-matrices in evaluation 
representations. The analytic continuation in all the parameters of the expression ( 4.11| ) 
can be given in terms of gamma and hypergeometric functions: 

-fc.i _ / 1 \k+n _ n sin P ~ ^] T ( l + !) 



= (-l) fe+n 7rL> 



sin[4vr] sin[(/s + I - £ 2 - n)ir] T(l -£ 2 + l)T(n + 1) 



r(n + 1 - £i)r l + - n - Su ) T ( 1 - - 5u 



r|Hi-^-fc + i r ^-fa 



x (4.14) 



4F3 ( —k, —n, 5u + 1 -, — — 5u; 1 — £1, £ 2 — k — I, I — n + 1; 1 I , 



where one has defined ±F$(x, y, z, t; r, v, w; r) = 4^3(0;, y, z, t; r, v, w; T)/[T(r)T(v)T(w)]. 

Moreover, we can easily see that we are in a special situation, since the parame- 
ters entering the hypergeometric function 4^3(01, 02, 03, 04; 61, 6 2 , 63; 1) satisfy Ylt=i a i ~ 
Sj=i bj = —1- When this happens, the hypergeometric function reduces to a 6j'-symbol, 
according to the following formula (see for example [56]): 



P , , , , n (-i) fel+1 r (62) r (63) y/r (1 - jjg r (1 - gj r (i - jg 

4F3 (ai, o 2 , 03, 04; &i,6 2 , 63; 1) = — , 

T (1 - 61) yT (6 2 - ai) T (6 2 - a 2 ) 

(1 - q 4 ) r (ai - 61 + 1) T (a 2 - bi + 1) r (a 3 - h + 1) r (a 4 + 

/ — 

A/r (6 2 - 03) r (6 2 - 04) r (6 3 - 01) r (6 3 - o 2 ) r (& 3 - o 3 ) r (6 3 - a 4 ) 



\ (-01 - a 4 + 63 - 1) i (-01 -a 3 + 6 2 -l) 2 (01 + a 2 
1 (-02 - 03 + 63 - 1) \ (-a 2 -04 + 62-1) 5 (03 + a 4 - 

By identifying the parameters we see that the relevant 6j-symbol 

J h h h 1 
1 J4 J5 J6 I 




(4.15) 



(4.16) 
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has coefficients 



3l = 2\ 2 

* = ^-,-i-* 

is = j (*i - 2 - fc - «) > 

J5 = 2 ( — 2 l-«-i + n + die 

k = \{i2-l)- (4.17) 

For generic values of Su, the 6j-symbol is understood in the sense of analytic continuation. 
However, one can prove that, for values of Su corresponding to the physical poles (see also 
the discussion below), the entries of the 6j-symbol are indeed half-integer, as one may 
expect from the fusion rules of su(2) representations. 

In the special case I = (a similar argument would hold for k = 0), we can go back to 
expression (|4lll ), and see that it can be casted in the following form: 



4.2 Poles 



r(fc + l)r(l + n - £ 2 )T (l - ^ - SujT (k + ^ - ^ - n - 5uj 

r(i - e 2 )r{k -n + i)r(n + 1) r (k - - + 1) r - 5u 



Next, we will analyze the pole structure of formula ( 4.11 ). For simplicity, we will restrict to 
the special case of expression ( 4. 13| ) in this section. The general case will then be analyzed 
in Appendix |B[ 

For the remainder of this section, we rescale the coupling constant and the spectral 
parameters according to g — > gV2, u — * r^u, in order to adapt our conventions to those 
of [48]. Formula ( f4.13| ) becomes 

-k,o _ rk\UU(h-plU k ^n-^[2 q -h + i. - 



« = d I — ^-7^ ; J v \ TV, • ( 4 - 18 ) 



We can recognize, following [48], the presence of potential poles of this formula at bound- 
state rank 2p — 1\ — £ 2 , which we now want to study. To begin with, in order for these poles 
to be 'physical' (namely, with a positive bound-state rank), one must have 2p > l\ + l 2 
(inclusion of the "equal" case will not affect the result). Since p is at most as large as 
k, and, from the definition ( 4.10| ), l\ > k + 1, we conclude that one needs i\ > l 2 + 2 



and 2£\ — 2 > 2k > l\ + t 2 . If this holds, then the physical poles occur for the values of 
k > p > (li + i 2 )/2. There are two possibilities, which we analyze in what follows. 



- 18 - 



The numerator does not identically vanish, and does not have zeroes for any physical 
values of the spectral parameters. This situation would leave physical poles uncan- 
celled in the denominator. This happens if £2 > n + 1 (by looking at the product 
]lp'=i(^2 - p')), and if h - £ 2 > 2{k - n), with k > n from ( pLlQl) . Combining the 
two conditions, we obtain k < n+ [(£1 —£2)/ 2] < — 1 + \{£\ +£2)/ 2], which contradicts 
the original condition for physical poles given above. Therefore, this situation cannot 
occur. 

The numerator does not identically vanish, but it has zeroes at some physical values 
of the spectral parameters. The latter can potentially cancel some of the physical 
poles, and we want to see whether few poles will be left uncancelled, or all of them 
will be neutralized. In order to have zeroes at physical values, we need £\ — £2 < 
2(k — n) — 1. This means that the interesting zeroes will occur at positive bound-state 
ranks running from up to — k + [{£\ + Z%)/2\, while the interesting poles occur at 
bound-state ranks running from up to n - k + l + [{£\ — £2) /2]. But we see that the 
number of these zeroes is always larger than the number of physical poles. In fact, if 
one subtracts the two numbers, one gets £2 — n — 1, which is bigger than (or at least 
equal to) zero in order for the amplitude not to identically vanish (once again, by 
looking at the product 1X^=1 (^2 ~p'))- Therefore, also this second situation cannot 
occur, and we conclude that we cannot have physical poles in formula ( 4.15| ). 



The result of this section and of Appendix |B] is consistent with the fact that in the su(l|l) 
sector corresponding to Case I, one does not expect any physical bound-state poles [16, 48]. 
The factor of D in fact cancels the s-channel pole at bound-state rank £\ + £2 coming from 
the overall scalar factor, and no physical poles are left in this amplitude. 

4.3 The classical limit 

In this section, we want to take the classical limit [57] of the S-matrix for scattering 
of arbitrary states belonging to Case I, and show that it is reproduced by the universal 
formula given in [52]. For general transition amplitudes, though only in the situation where 
one has at most two bound-state components, this test has already been performed in [47]. 



When looking at formula (4.11), one can see that, besides expanding the factor D, 
one needs to expand the remaining expression, depending only on the difference 5u of the 
spectral parameters, for large values of Su. In the classical (near BMN) limit, in fact, the 



coupling g goes to infinity 5 . With our conventions ( 2.26j ), the spectral parameter u grows 



linearly with g in this regime, while both parameters x , expressed as 



tend to their common classical value x [58]. 

The relevant terms to the classical limit of ( fi.ll ) are given by the following expansion: 



3 The momentum goes then to zero as p ~ l/g. 



- 19 - 



f 1 - h - p)) m=i(h - i) nil- n ^ - o 

~ (1 + L> d )- r ^ ; X 




A; + i — m — p+2 
2 



(4.20) 



where D c ; denotes the first order in 1/g of D. Here, we have used the fact that the binomials 
enforce I > n — m, in order to obtain the power of Su 2m ~ k ~ n . Let us start by considering 
non-diagonal amplitudes, namely, n different from k (cfr. ( [4.10| )). In order to do that, let 
us first reduce the above formula for the case n > k. In this case, the leading piece in the 
above expression is given by the term in the sum with m = k (the binomials are in this 
case non-zero, since, from (4.10), one has I > n — k). The amplitude tends to 



As one can see, in the non-diagonal case only one of these amplitudes actually contributes to 
the classical limit (corresponding to the order 1/g of the scattering matrix). Namely, only 
the transition from a state characterized by quantum number k to one with corresponding 
quantum number n = k + 1 has the right order, the other ones being suppressed by higher 
powers of 8u. In this situation, the classical amplitudes reads 

We checked that this is exactly the value one gets from applying the universal formula of 
[52] to this amplitude. For convenience of the reader, we report here below their classical 
r-matrix. 

Ti2-S®M-1HI®S £®H BIOS i, , 
ri2 = — — 77-^ r-^ — + — — + -K - u± )M ® M, 4.23 

l(U\ — U2) IU2 lUl 2 



with 



and 



T 12 = 2 ( R a <g> - L b a ® L a 6 + G a Q ® Q a a - Q a a ® <G a a ) (4.24) 



~ 9 a — ) . (4.25) 



2 ad + be \ a dw a a d9 

In ( 4.23j ), we understand all generators (taken at their classical value) as differential 
operators, acting on the appropriate monomials corresponding to Case I states. We then 



- 20 - 



compare the result with the expression we have found above for the first order in 1/g of 
the complete S-matrix. 

Next, let us consider k > n. In this case the binomials force the leading piece in the 
sum to be the one with m = n. This reads (quite symmetrically w.r.t the previous case) 

ark,. 1 Y[U{Zl-i)Y[tl- n ^-i) ( k \ 

n ^ k ~ n iiid-ionU^-p) \*- n /£± ~ =fts " 

Analogously, only one of the non-diagonal terms has the right falloff to be able to contribute 
to the classical r-matrix, namely the amplitude for quantum numbers k to n = k — 1. The 
contribution is given by 

~ ^k(£ 2 -l-l), (4.27) 

and we also checked it against the classical proposal of [52]. 

The diagonal part, for n = k, is slightly more complicated. The leading term can be 
obtained by specializing to k = n either of the two formulas ( [4.21 ) or ( 4.26 ), and is easily 



seen to be equal to 1, as expected. The quantum R-matrix goes in fact to the identity in 
the strict classical limit. The next to leading term of order 1/Su contributes to the classical 



r-matrix, and can be straightforwardly obtained from (4.20) as 
1 



_p=l p=l p=l— k 



This expression can be explicitly evaluated, and, after supplementing it with the suitable 
overall scalar factor (cfr. [47]), we have checked that it precisely corresponds to the result 
coming from the universal formula of [52]. 

As a curiosity, we have checked that the order 0(l/g 2 ) of the Case I amplitude is 
completely reproduced by half the square of the classical r-matrix. The departure from a 
simple exponential series seems to reveal itself starting from the next order 0(l/g 3 ). We 
plan to come back to this issue in the future, in the light of possible consequences for the 
abstract form of the universal R-matrix (cf. for instance [59, 60]). 

5. The S-matrix for Case II 

In this section, we will use the S-matrix derived for Case I to find the S-matrix for Case II. 

As explained in the previous sections, AQg, AG| and their Yangian partners map Case 
II states onto Case I states. We introduce the Case II S-matrix in the following way 

S M? = EE ^ j \m,N-rn)f, (5.1) 

j=l m=0 

where again N = k + I. This means that the coefficients & m if actually correspond to the 
S-matrix restricted to the following spaces 

: Kl — O-n- (5-2) 
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Generically, both spaces are 4 dimensional, and correspond to the coefficients of a 

4x4 matrix. One might wonder what happens for n = 0,N since e.g. |0,iV)^, strictly 
speaking, does not exist. However, it turns out that these states are always multiplied by 
0. Hence, the 4x4 matrix actually contains the non-generic case in which n = 0, N. This 
will be explained later on in Section and we will continue with deriving the generic 4x4 
matrix. 

By now considering the action of A op Q^, we can relate the Case II S-matrix to ( 4.1 1] ) . 
It is easily checked that 



l\k,l)Y = Q j (k,l)\k,Q I , 



with 



(5.3) 
(5.4) 



Qi(k,l) = a 2 (l-£ 2 ), Q 2 (k,l) = ai(h-k), 
Q 3 (k,l) = b 2 , Q A {k,l) = -b 1 . 

Similar expressions are of course obtained for A op (Q>3, A op G|, AG|. We can now apply our 
general strategy in the following fashion: 

4 k+l 



(n,N- n\ A op Q^§ \k,l)f = J2Y1 N ~ n ' A ° PQ 3 ' m ' N ~ m ^ 

j=l m=0 
4 k+l 



Yl Yl N - m ) l ( n > N - n i m ' N - m ) 1 

j=l m=0 



n) 



3=1 

On the other hand, we can use the symmetry properties of the S-matrix to obtain 

l (n, N-n\ A op QlS \k, l)f = l (n, N — n\ SAQ3 \k, l)f 

= Qi{k,i) l {n,N -n\ § \k,l) n 



(5.5) 



N 



Qi{k, I) &m N - n K N - m ) 1 
Qi (k, V) S£ n ' ■ 



(5.6) 



Clearly, this gives us four linear equations relating the two S-matrices. A similar compu- 
tation can be worked out using A op G 2 , giving four additional equations. We can cast the 
above formulae in a convenient matrix form: 



ai(N— n— £2) az(h—n) 64 —63 
C4(N— n— £2) C3(£i—n) di — dz 



k.l 



k.l I a 2(l-k) ai(£i—k) 62 -bi 
c 2 (i-fe) ci(£i-k) d 2 -di 



(5.7) 



with 



k.l 



1 ay k ,l;l 


^n;2 


<v,k,l;l 


• y n;4 


^k.l-,2 
- y n;l 






• y n;4 


a^k,l;3 

- y n;l 




<3 !;3 


twk,l;3 
• y n;4 


V»&* 4 


■ y n;2 


ayk,l;4 
,y n;3 





(5.8) 
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Written in this way, the relation to (|3.l| ) becomes apparent. However, it is clear from the 
above matrix equation that, in order to fully determine &n' 1 (and therefore the full Case 
II S-matrix ), one needs additional equations. 

These equations can be obtained via the Yangian generators. Consider the following 
operators: 



AQ = A(Qi) + 



AG = A(G|) + 



2AL^A( 



h + i 2 - 2(N + 1 + 5u) 
2AL^A(Gf) 



+ 



i- l 2 + 2(N -2n + Ul + u 2 ) 
2(h + h) - 4(iV + 1 + 5u) 



h + 2(N - 2n + m + u 2 ) 



1 + £ 2 -2(N + l + 5u) 2{£ 1 + t 2 ) - 4(N + 1 + 5u) 



ALiA(Q|), 

(5.9) 
AL*A(G?). 



These operators are chosen in such a way that only states of the form \n,N 



n > 



arc 



mapped to |n, N — n) 

this is important in (5.5) in order to be able to pull out the matrix 
generically maps 



k.i 



AQ : V; 

or, more precisely, we can write 
AQ|M>" = 



When we follow the same derivation as in the above, we see that 

In fact, AQ 

VU,l+V (5-10) 

(5.11) 



ii 

k,l 



k+i,i-i 



V, 



k,l 



Qi(k,l)\k,l) 1 + Qf(k + l,l- l)\k + 1,1- l) 1 + Qi(k -1,1 + l)\k -1,1 + 1) 
This clearly means that, if one follows Q5.5|), one obtains 



l (n,N-n\ A op QS \k,l)f 



(5.12) 



4 



^.f d Qf(n,N-n) + ^ 



k,l;j Aop+ 



(n, N-n) + (n, N — n) 



However, the specific choice we made for AQ means that Q° p ' + (n, N — n) = Q° p '~(n, N — 
n) = 0. In other words, we can again pull out the matrix factor ty*' 1 on the left hand side 
of the final equation. Since this is specifically tuned to work for the opposite coproducts, 
the right hand side of the equation will not have this property, and will contribute 
there. This is exemplified in Figure [|. 

For compactness, let us define M = N — 2n (= N — n — n). By combining all the 
equations one is lead to the following matrix equation: 

a 3 \ (a 2 ai \ 



C4 C3 
«4 03 

\o c 4 C3 y 



k.i 



C2 ci 
«2 «i 

\o c 2 Cl y 



B + ^ +1 ' 1 - 1 + B- + B%^ 1 \ 



where the matrix on the left hand side is given by 

/ N-n-h 
n-t\ 



A 



>34 
^34 
1 

■S43 



"^43 

\ (JV-n-€ 2 )(<i^ 43 ) (n-4)(5«+M) m-Su^+^-s^^ 



1 

^34 
>43 

^34 _ 

j+Af+l 1 -l 2 - g 34^34 
^43 
(Su+M+lj) ^43 



(5.13) 



34 
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A op Q, A° P G 




Figure 2: Action of A op Q, A op G and AQ, AG. 
states (on the right). 



AQ, AG 




They map Case II states (on the left) to Case I 



and the matrices on the right hand side by 



B+ 



2{l x - k — l)c 



l—n 



-N 



/o o o o\ 



•±Z12 

V 00 <Al °J 



/o 



2{l 2 -l-l)c 



n—l 



-N 








^— 



B 



l- 



&12 

1 

&21 
(N-8u-l 2 ).J? 12 

&12 



_JL_ \ 

^21 » 

■^21 
^12 





i\-k 
(l-l 2 )(N-8u) {tx-k)t 2 J?]_ 2 
\(h-l)(h^2i) {h-k)(5u-N) '-"Wi2+'a (fa-iv+W 21 / 

/ 



JV-iSu-l]^ -^2-^12^1 
^12 

1 12 



-JV 











\ 



o 

i(l+n+fc-£i) (/-^ 2 )(l+n+fc-£ 1 ).y 12 -(^i-fc)(l+Af-n+i-£ 2 ) 

fe(l+iV-n+«-fe) -(/-£ 2 )(l+n+fc-<?i) ^i-fe)(l+Ar-n+«-fe)^2l/ 



where we defined 



is ij — Ctidj bjCi. 



(5.14) 



These coefficients satisfy the following identity 



(5.15) 



Notice the similarities between the matrices A, B and B + , B~ . From this, it is now straight- 
forward to extract ' l by simple linear algebra. For completeness, we will explicitly report 
this inverse matrix. Define 



Da = det A 



6)2 



(n-4)(iV-rt-4) 



-45U 2 + (£l - l 2 f + 44^34^43 
-40+^+4^4^34^43, 



(5.16) 
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then 



A" 1 = — 
da 



I 



\ 



f -f+[M+ 



N-n-e 












1 

n-e-i 



X 



\ 


<S 4 3 
^43/ 

[tf +2^1^34^43] ^34(^+[M+^2]-C+) ^+^34^43 J+^34 \ 



(5.17) 



.^[M+^]^43 



[cr+^2^34^43]-^ 



€2 [M+^]^34 



Cfj*43 
2<?1^43 

— cT 



h +^2 ^34^43 



-2fe.-^34 / 



and 



k.l 



A' 



Sl32 #31 

B 3A ^34 

^42 ^41 

^43 ^43 



\ 













Stv. 







V 



_^ £31 

S34 i2 34 
^42 ^41 / 
^43 ^43 / 



{ t r n fe+1 ' z - 1 s+ + ^; A; - 1 ' / + 1 i?- 



+ .r^'fl} . (5.18) 



Note that the final result for & n ' purely depends on the spectral parameters through their 
difference 6u, and the representation parameters only appear in the combinations £?ij,7iij 
(modulo perhaps an overall factor), the rest being taken care of by combinatorial factors 
involving the integer bound-state components. 



6. Complete Solution of Case III 

We will perform here a similar construction as done in the previous section, in order to 
solve Case III in terms of Case II. Let us first set few additional notations. We introduce 
the S-matrix at this level in the following way: 



S\k,l) 



111 



k+l 6 



2.2. 

m=0 j=l 



J 



(6.1) 



It is clear that one can repeat a very similar derivation as performed in ( |5.5|) and fl5.6j), 
where, instead of S£ , one has to think of having <3f (and indices running over the appropriate 
domains). This time, moreover, one considers the action of AQ3, AG|. The result is now 
the following matrix equations: 



/ (n-£i)a 3 
(N-n-£ 2 )a4 64 

(n-£i)a 3 




h 





b 4 

b s (N-n~l 2 )a 4 






(I— ^2)02 ?>2 

(fc-4)ai 



(N-n— £2)14, 64 

61 b 2 



61 (i~e 2 )a 2 











(i~e 2 )a 2 b 2 



\ 

-63 



(n-£ 1 )a 3 J 

\ 

-61 





k.l 



(6.2) 



(6.3) 
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and 



/ (n-£i)c 3 d 3 

(N~n-£ 2 )c 4 d 4 
(n-£i)c 3 



where 



d 4 


d 3 {N-n-e 2 )d 4 





-d 3 




(N~n-e 2 )c 4 d 4 







J, 



k.l 



y 



k.l 



di 




dz 






(l-i 2 )c 2 d 2 

o (fc-^i)ci o di (i~e 2 )c 2 

(Z-fe)c 2 d 2 



(6.4) 



o 
-di 
o 

(fe-£l)ci/ 



^ i 



k.l 





"^n;2 


,1; 1 
-*n;3 






tiwk ,1 ;2 
"^n;2 


,Z;2 


^fc,Z;2 


■^■n-,1 


,Z;3 
"^n:2 


opk ,1 ;3 
"^n;3 


<5¥>/v,Z;3 


^n;l 


"^n;2 


-*n;3 


™fc,Z;4 
-^n;4 




^yfc ,Z;5 
"^n;2 


^jfe ,Z;5 
-^n;3 


^fc,Z;5 
-^n;4 




-*n;2 


,Z;6 
-^n;3 


^fc,Z;6 
"^n;4 



fc,Z;l 
n;5 

fc,i;2 
n;5 

fc,i;3 
n;5 

fc,i;4 
n;5 

fe,i;5 
n;5 

fc,i;6 
n;5 



M;i\ 

;6 \ 

c,Z;2 



(6.5) 



Again the relation with (3.1) is apparent. 

However, it is readily checked that these equations are not independent. Hence, once 
again one needs additional equations, as in the previous section it was required in order 
to compute <3/ . In that case, they were provided by Yangian generators. In this case we 
are more fortunate and do not need the Yangian, since one can consider the action of AQ4 
and AG?. It is easy to check that, by repeating the above procedure using these additional 
symmetries, one arrives this time at the following matrix equations: 



( na 3 
(N—n)a 4 64 
na 3 



63 












63 




V 







(N—n)a 4 64 — na 3 



-64 


(iV— 71)04 




k.l 



k.l 







61 
la 2 b 2 
fcoi 61 
\0 la 2 b 2 




61 


-kai 



-b 2 \ 


/«2 
/ 



(6.6) 



and 



nc 3 
(N—n)c 4 d 4 
nc 3 




d 3 








d 3 



(N-n)c 4 d 4 

where we have defined 





d 3 

-nc 3 



—d 4 


(N-n)c 4 




^ 'feci 



6£ 



k.l 



k.l 







lc 2 d 2 
fcC! 



d 1 

d 1 




di 








lc 2 d 2 —kc\ 



-d 2 \ 


lc 2 

/ 



(6.7) 



rtfk,l-\\\ frisk 



hi 



'n-l;l 



^n-l;3 ^n-VA \ 



fi>/K— 1 ,/;2 [j/skj, — 1;2 (jj/k— 1,Z;2 (j//k,l— 1;2 
°*n;l ^n;2 ^n;3 ^n;4 



(ji/k— 1,Z;3 (j)/k,l— 1;3 (jj/k l,i;3 
^n-l;! ^n-l;2 ^n-l:3 



1 ^fc-l,/^ ^fc,Z-l;4 ^fc-l,Z;4 ^fc,Z-l;4 



(6.8) 



n;3 



n:4 
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Combining all of the above equations is sufficient in order to solve for 3? . To be more 



precise, one can write the equation for 3fn' 1 in the following way: 



/(n- 







^43 





1 


o \ 







1 











-^43 












^43 















~J*43 








1 







-1 








-J*43 





V 





-71.S43 





-^43 





o / 



yk,l 



(6.9) 



k.i 



( (k 







^32 





S*2 





\ 




-£ 2 )^42 


J*42 











— ^32 









(fc-£l)^32 





^32 


(«-fe)^24 















(/-£ 2 )^42 


J*42 





{k-tx)B32 






kg 2 3 





-J%2 








,J?42 






1^24 


-^42 








-^32 












k^23 










^24 




V 








/<@24 


-y±2 


fc^31 





/ 



with 





•*n;2 


ayk,l;l 
,y n;3 


ayk,l;l 

' y n;4 


















\ 


rjyk,l;2 
' y n;l 


ayk,l;2 
■ y n;2 


,wfc,i;2 
"n;3 


^fc,!;2 
■^4 




















^fe,(;3 
• y n;l 


awk,l\3 
■ y n;2 


n;3 


Oi/k,l;2 
































<wfc-l,i;l 
-J n-l;l 


ayk,l- 

' y n-l 


1:1 

2 


J n- 


-l,i;l 
-1;3 


0>/k,l-l;l 
■ y n-l;4 
















<wfc-l,i;2 


ayk,l- 

■ 9 n;2 


-1;2 


at/k 
n; 


-l,i;2 
3 


iwk,i-ia 

•^n;4 




V 











<wfc-l,i;3 
' y n-l;l 


ayk,l- 


-1;3 
2 


W k 
•"n- 


-l,i;3 
-1;3 


<wfc,!-l;3 
■ y n-l;4 


/ 



(6.10) 



The explicit matrix inversion gives 
/ 1 



avk,l 








£l^34 ^1^34^43 



n n-^i 

^1^43 









1 



1 



1 







1 

' -^43 



£l-S 3 4 
1 

^34^4 2 3 
n 

1 

£1^43^43 
1 

£1^34^43 
/ (k-h)J3:n 

(1-^41 




fci2is 
U2i4 





£l<S 3 4 £1^34.-^43 



n— £1 n 
hJ?43 J^Jf 3 







1 

' -^43 



£l^34 
1 

£1^34^3 
n-ij 

£l^43 
1 

£1^43^43 
_ 1 

£1^34^43 



k,l 



(6.11) 






^31 





S41 





^41 











-^31 







.^31 


(«-£ 2 )^14 








(/-£ 2 )^41 


S41 





(fc-^l)^Sl 





-.^31 








J*41 


-Sai 








-S31 





kSi3 





-S31 





J.3l4 





/J2l4 




kJ2 3 i 






It is now straightforward to do the matrix multiplication. This solves the final case. Once 
again, the dependence of the entries solely on the difference of the spectral parameters, 
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and on the characteristic combinations of representation labels already observed in Case 
II, is a noticeable feature of the result. 



7. Reduction and Comparison 

Let us now compare our formulae with the known S-matrices. Here, one runs into potential 
difficulties. The formulae from the previous sections were derived for generic bound states, 
and one might wonder whether there could be obstructions for small bound states. A 
first problem arises when n is comparable to £\ , £2- A second problem is encountered for 
n = 0, n = k+l, since the basis of two-particle states in these two cases is lower-dimensional. 
One can wonder whether our formulae 

k+l 

§\k,l) 1 = J2^'n l \n,N-n) 1 (7.1) 

n=0 
k+l 4 

S\k,l)f = J2J2^t' d \n,N-n)f (7.2) 

n=0 j=l 
k+l 6 

SIM)? 1 = EE^i"- JV - ft ?- ( 7 - 3 ) 

n=0 j=l 



with N = k + l and 3? given by (fOJD and ( jsTTH ), remain valid also for these particular 
values. 

It turns out that this is indeed the case. Let us deal with the first problem. One 
can see from ( 4,11| ) that, when n > £\, precisely the unwanted S-matrix elements vanish, 



basically thanks to the vanishing of the correspondent coefficients 3t>n' 1 ■ 

Concerning the second potential problem, we notice that the issue arises only for Case 
II and III states. In Case II, the corresponding sum on the right hand side of ( [7.2] ) contains 
terms like 

r k,lA 



%f A %N) l l (7.4) 

But, as seen from ( |2.14D , |0, N) 1 ^ is not well-defined (actually it is not part of our bound 
state representation). Hence, the S-matrix transition amplitudes toward these states, 
$Q fe /' 4 , should vanish identically. We verified that this indeed turns out to be the case, 
which means that these states completely decouple. 
More specifically, from (|4.11 ) it can be shown that 



ay-k+1,1-1 _ h I ay-k,l , 7 fl 

Su - - I 



= — l——%- Q k > 1 . (7.6) 



This means that in ( 5.18 ) one can pull out a factor . The remaining matrix part is 

fk,i 
0;i 



straightforwardly seen to have zeroes for the states corresponding to the amplitudes <3C!f' 1 1 
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for all i = 1,2,3,4 as indeed should be the case. In other words, one can unambiguously 
write 



k+l 4 



n) 



(7.7) 



n=0 j=l 



where ty^.j^ is given by the complete 4x4 matrix from (5.18). The same should be true 



n;i 



for Case III states. 

One can now compare our coefficients against the known S-matrices. Complete agree- 
ment is found with E> AA , S AB , S BB from [6, 15,43]. We also checked several coefficients of 
the S-matrix S 1 ^ from [31] which also agree with our findings. 
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A. Yangians and Coproducts 

The double Yangian DY(q) of a (simple) Lie algebra q is a deformation of the universal 
enveloping algebra U(q[u,u^ 1 ]) of the loop algebra g[u, u^ 1 ]. The Yangian is generated by 
level n generators S A , n£Z that satisfy the commutation relations 

KX\ = ^ B jg +n + o{h), (a.i) 

where F AB are the structure constants of g. The level-0 generators § A span the Lie-algebra. 
The su(2|2) Yangian can be supplied with a coproduct in the following way [44-46] 

h n ^ 

A(J ? f ) =J>1+ UM ® Jtf + - £ FBct-i- m U [C] + 0(h 2 ), 

m=Q 

A{U)=U®U, (A.2) 

where U, the 'braiding factor', equals e*2 , and h = 1/g. 

An important representation of the Yangian is the evaluation representation. This 
representation consists of states |tt), with action l A \u) = u n S A \u). In this representation 
the coproduct structure is fixed in terms of the coproducts of Jo, Ji- For the remainder of 
this paper we will work in this representation, and identify Ji = J = nj for the su(2|2) 
Yangian. The spectral parameter u depends on x^, as one can see from formula ( |2.26| ). 

The S-matrix is a map between the following representations: 

§: nM,jn®n,(P2A)^viApi^)®nAP2,j pi i (a.3) 
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where 'fe i (pi,e 2 ' 1 ^) is the £j-bound state representation with parameters a,i,bi,Ci,di with 
the choice of £ = e 2 ^. This specific choice removes the braiding factor IA from appearing 
explicitly in the formulas [15]. 

The bound state S-matrices are now fixed, up to an overall phase, by requiring invari- 
ance under the coproducts of the (Yangian) symmetry generators 

S A(JJ A ) = A op (J A ) S, 

S A(J A ) = A op (J A ) S, (A.4) 

where A op = PA, with P the graded permutation. 

For completeness and future reference, we give the explicit formulas for the coproducts 
and for the parameters aj,6j,Cj,dj. First, the su(2|2) operators: 

A(jtf) = j£ +j[£ . (A.5) 

The coproducts of the Yangian generators are then given by [44] 

A(L a 6 ) = £,£ + l 2 % + - \^ 2 % ~ l G ^My ~ \<*& G & 

A(R^) = R£ + ±$ - \r^R 2 % + + \g,%Q 2 % + \^G 2 % 

-ls$Gj c Q& - ^Q&G 2 7 C , (A.6) 
A(Q« /3 ) = 0i5, + - ^Q 2 ? 7 + ^Qi^J " ^Kc Q2^J + 

A(G Q b ) = 6$ + G 2 f b + \l^G& - \g x %1. 2 % + ^ G 2 J - ^JR^ 
+\m 1 G 2 % - Ig$K 2 - \e bc e^C\q 2 ^ + \e bc e^Q^C\ , 

and for the central charges 

A(H) = Mi + H 2 + CiC 2 - CjC 2 , 

A(C) = Ci + C 2 + ^HiC 2 - ^CiM 2 , (A.7) 
A(C f ) = Cj + C 2 + ^HiC^ - iclM 2 . 

The product is ordered, e.g. QiQ 2 means first applying Q 2 , then Qi (as differential oper- 
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ators). Finally, the coefficients used in A are given by: 



*P2 / 9 Vi 



Cl 



e 4 e 2 



i • 



°2 = y^m, 

02 ~ V 2* 2 x + 

■ P2 I 

r]2 = e l -Jix2 



h 

di 



b 2 
d 2 



-te 



W2 



I -1. J_ ( £L _ 1 



2*i m \xf 



1 



'_S_J_ i Ed. _ i 

2*2 IB Uj 



I 9 x 2 ( x 2_ 
2*2 M?2 



(A. 



The coefficients in A op are given by: 



03 



C3 = 
op 



Mi T 'l ' 

op 

9 Vi 



e 4 



2X 



1 • 



d 3 



2*1^7 



04 



op 
V 2 



_9_°V 

2l 2 '2 

-ipi 



op 

_9_!hL 

2*2 x+j 



e 4 e 2 



b 4 
C?4 



-?.e 



9 1 
2*2 5g" 



1 , 



2*2 ig* 



The non-trivial braiding factors are all hidden in the parameters of the four representations 
involved. 



B. Poles of the general Case I Amplitude 



Let us now turn to formula ( [4.1l| ), and outline a proof that it also does not have poles 
at physical values of the momenta. In order to do this, we will make use of its rewriting 
in terms of the 6j-symbol according to formulas ( 4.15| ), ( 4.17| ). The advantage is that the 
6j-symbol itself is analytic, and the singularities are essentially factored out in terms of 
gamma functions. 

First, by looking at ( 4.11|) , and remembering the discussion of the case I = 0, we see 
that the source of possible physical poles is the denominator 1 / Hpi| {$ u + ^j^ 2 — p) ■ In 
order to have poles at a positive (nonzero) bound-state rank we need (see Section \l.2\ ) 5u = 
p-[(£ 1 +£ 2 )/2],p€ ((h+h)/2, k+l], k+l > (h+£ 2 )/2. This denominator is recognized as 
the ratio of gamma functions T ( / + 



*l-*2 



n 



5u 



/r(k + i 



*i+* 



Su + 1 ) in flP3 ) 



From ( 4.11 ), it is clear that the only thing that can happen is the rest of the formula 
cancelling some of these poles with zeros. Let us first analyze just the contribution coming 
from the hyper geometric function, rewritten as in ( [4.15 ), After identifying the values of 
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the parameters, namely 



a\ = —k, a 2 = —n, 

a 3 = 1 1- du, 

4 - ii , 
a 4 = Su, 

61 = 1-4, b 2 = -k-i + e 2 , 

6 3 = l + /-n, (B.l) 
one can see that four of the gamma functions in the numerator and four in the denom- 



inator of ( 4.15 ) bear dependence on Su. At the values of Su corresponding to physi- 
cal poles, two of these gamma functions in the denominator become T {—k — I — p) and 
r (— 1 — k — l + 1\ + £ 2 — p)- The first one has always negative argument, since p < k + I. 
The second one as well, since -l-k-l+£i+£ 2 -p < -l-k-l+[(£ 1 +£ 2 )/2] < -1-2&-2Z < 
0, where we have used the physical conditions on p and k + I reported above. These two 
poles under the square root combine to give a simple pole in the denominator, namely a 
zero, that cancels the physical pole coming from the ratio of gamma functions discussed 
above. However, this does not terminate our analysis, since we still have to make sure no 
other physical poles are generated from all the other gamma functions still present in the 
formula. Let us spell those singularities out. 

Two of the four gamma functions depending on Su in the numerator are always regular 
at the physical value of momenta. They reduce in fact to r (1 + p) and V (£± + £ 2 — p). 
The first one has manifestly positive nonzero argument, and so does the second one, since 
£\ + £2 — P > £\ + £2 — k — I > 2 using conditions on p and £\ > k + 1, £ 2 > I + 1 from 



(4.10). Poles could only arise from the remaining two gamma functions depending on 



Su in the numerator, which at the physical values of momenta become T (1 — £ 2 — p) and 
r {£\ — p). These two poles are however cancelled by a zero coming from T - Su) m 



the denominator in ( 4.14 ). The latter reduces to T (£1 — p) at physical momenta. Factoring 
it into two identical square roots, we see that one of them exactly cancels one contribution 
from the hypergeometric function, and the other one cancels the singularities of the other. 
In fact, when r (1 — £ 2 — p) has poles, namely for 1 — £ 2 + p < 0, this implies £ 2 — 1 > 
p > (£1+ £ 2 )/2, i.e. {£1 - £ 2 )/2 < -1. But this means that l x - p < (£1 - £ 2 )/2 < -1, 
and the two poles cancel. In order to conclude the argument, we still need to analyze 
potential poles coming from the gamma function r ( I + ^"^ 2 — n — Su) in the numerator 



of (4.14). But when this gamma has a pole, namely for I + £\ — n — p < 0, one of the 
four gammas depending on u\ — u 2 in the denominator of formula ( f4.15| ) (and different 
from the two already considered), becomes r (I + £\ — n — p), which has also a pole. So 
the latter reduces the upper pole to a square root of it. However, the original formula does 
not have any branch cut, therefore it is not possible to have one square root singularity 
left uncancelled. Since one can check that all the other parts of the formula, excluding the 
6j-symbol, are neutralized, one concludes that the 6j-symbol must have a zero in this case. 
We checked that this is the case for few examples, where one can see that the triangularity 
condition for the 6j-symbol is violated. 
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